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1.  Introduction  and  formulation 

It  is  well-known  that  there  can  be  a  significant  difference  in  behavior  in  flight  between 
liquid-filled  and  solid-filled  projectiles.  The  difference  is  caused  by  the  motion  of  the  liquid 
inside  the  spinning  projectile.  This  motion  causes  forces  to  act  on  the  projectile  which  can 
ultimately  cause  the  flight  of  the  projectile  to  be  prematurely  terminated  by  instability. 
The  initial  motion  of  the  projectile  necessarily  causes  the  fluid  motion  in  the  cylinder  to  be 
time-dependent;  later  it  can  be  assumed  that  the  flow  is  steady.  It  is  the  latter  situation 
which  is  investigated  here. 

Our  concern  then  is  with  the  motion  of  an  incompressible  viscous  fluid  of  kinematic 
viscosity  i/ina  cylinder  of  length  2 Aa  and  radius  a.  The  cylinder  is  spinning  at  a  constant 
rate  about  its  axis,  which,  in  turn,  is  nutating  at  a  constant  rate  and  yaw  angle  about 
an  axis  directed  along  its  trajectory.  Throughout  this  calculation  a  will  be  used  as  an 
appropriate  length  scale  while  the  magnitude  of  a  typical  velocity  is  taken  to  be  (ft  + 
f  cos  Ao)a,  where  ft  is  the  magnitude  of  the  angular  velocity  of  the  projectile  as  observed 
in  the  aeroballistic  reference  frame,  f  is  the  dimensional  coning  frequency,  and  A'o  is  the 
coning  angle.  (In  the  aeroballistic  frame,  one  cartesian  coordinate  lies  along  the  cylinder 
axis  and  a  second  coordinate  lies  in  the  plane  of  the  cylinder  axis  and  trajectory;  see,  e.g., 
Reference  1).  K0  is  assumed  small  throughout  this  work,  so  that  it  is  possible  to  determine 
the  properties  of  the  forced  motion  using  a  simple  perturbation  expansion  in  powers  of  A'o . 
We  nondimensionalize  the  coning  frequency  as  follows: 


r  = 


(ft  -I-  t cos  A'o) 


(1.1) 


and  we  define  the  Reynolds  number  by 


Re  = 


(fi  +  r  cos  A'o)a2 


(1.2) 


By  use  of  the  assumption  already  introduced,  A'o  <  1,  cos  A'o  is  replaced  by  unity  in  the 
above  and  in  the  following.  The  typical  time,  length,  velocity,  and  pressure  scales  used  to 
nondimensionalize  the  flow  variables  are  (ft  +  f)-1,a,(ft  +  f)a,  and  pa2(ft  4-  f)2,  where  p 
is  the  density  of  the  fluid. 

In  a  previous  paper,  Hall  et  al.I2!  ,  hereafter  referred  to  as  HSG,  we  used  the  spatial 
eigenvalue  approach  to  determine  the  forced  motion  in  the  projectile  at  finite  values  of  the 
Reynolds  number.  Here  we  aim  to  extend  the  results  of  HSG  to  the  high  Reynolds  number 
limit  and  show  that  in  this  case  there  is  a  significant  simplification  of  the  problem.  We 
now  indicate  how  the  appropriate  partial  differential  system  to  be  solved  can  be  derived 
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from  the  Navier-Stokes  equations  of  motion.  Our  derivation  is  of  course  only  brief  so  the 
reader  is  referred  to  HSG  and  the  references  therein  for  a  more  complete  account. 

First  the  Navier-Stokes  equations  are  written  down  with  respect  to  an  inertial  cylin¬ 
drical  polar  co-ordinate  system  (r,  8,x)  with  corresponding  velocity  components  ( u,v,w ) 
and  pressure  p.  The  coning  motion  of  the  projectile  means  that  the  boundary  conditions 
must  be  applied  at  a  moving  surface;  this  is  the  main  drawback  in  the  use  of  the  inertial 
reference  frame  but  for  small  angles  a  Taylor  series  expansion  of  the  different  flow  quan¬ 
tities  about  the  mean  position  of  the  cylinder  produces  boundary  conditions  at  a  fixed 
surface.  Thus  it  is  assumed  that  the  motion  of  the  cylinder  is  proportional  to  exp  i{ft  —  8) 
where  t  denotes  nondimensional  time  and  /  is  ,  in  general,  the  nondimensional  complex 
frequency  of  the  projectile  motion.  We  shall  take  /  to  be  real  so  that 

f  =  r.  (1.3) 

Since  7v0,  the  coning  angle,  has  been  assumed  small  we  expect  only  small  depar¬ 
tures  from  solid  body  rotation  and  the  velocity  components  are  therefore  (— A’ou*,r  — 
KoV*,—Koiv*)  and  the  pressure  is  |r2  —  Kop*.  Here  a  star  denotes  a  perturbation  quan¬ 
tity. 

It  is  assumed  in  our  calculation  that  the  motion  of  the  cylinder  is  specified;  we  do  not 
allow  for  a  feedback  from  the  fluid  to  the  projectile.  Following  HSG  it  is  convenient  to 
define  complex  velocities  and  pressures,  denoted  by  underlining  and  given  by 

( u*,v*,w*,p *)  =  Real{(u,y_,w,p)  exp i{ft  -  0}}.  (1.4) 

The  resulting  nondimensional  Navier-Stokes  equations  are 

*(/  -  Ik  ~2v  =  -v  +  Re-1  [V2u  -  -|*  + 

Kf  ~  Ik  +  2 u  =  (i/r)p  +  Ae_1[V2p  -  ~v  -  ^u], 

r ri 

i(f  -  l)w  =  -p  +  J?e-1  [V2u>  - 
(ru)r  —  iv  +  rwz  =  0, 

where  V2  =  52  -f-  +  d\  and  subscripts  denote  partial  derivatives. 

We  now  choose  to  seek  separable  solutions  of  these  equations  which  take  the  form 

( U,v,w,p )  =  (U(r)  sin  Kx,  V(r)  sin  Kx,  W(r)  cos  I\x,  P(r)  sin  I\x),  (1.5) 
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(1.6a,  6,  c,  d) 


and  in  this  case  it  ff'Hows  that  ( U ,  V,  W,  P)  satisfies 

[i?e'  \V71  -  r"2)  -  iM\U  +  2(1  4-  (i/r2Re)]V  -  Pr  =  0 

iP 

(i?e_1(Vi  -  r"2)  -  iM]V  -  2(1  +  (i/r2Re)]U  +  —  =  0 

-  iM)W  —  I\P  =  0 
(rll)r  -  iV  -  KrW  =  0. 

where  M  —  f  —  1 
and  Vi  =  drr  +  ^dr  —  [^r  4-  A”2] . 

These  equations  must  be  solved  subject  to  an  inhomogeneous  condition  at  r  =  1  and  an 
appropriate  condition  at  r  =  0.  We  obtain 

U-iV  =  W  =  P  =  0,  r  =  0  (1.7) 


at  the  centre  of  the  cylinder.  For  a  discussion  of  the  origin  of  (1.7)  the  reader  is  referred 
to  Batchelor  and  Gill^.  The  ordinary  differential  system  specified  by  (1.6), (1.7)  and  the 
the  no-slip  condition  at  the  cylinder  constitutes  an  eigenvalue  problem 

K  =  K(Re,f).  (1.8) 


At  finite  values  of  Re  there  are  no  real  solutions  of  that  eigenvalue  problem  so  that 
solid  body  rotation  is  stable  to  the  first  non-axisymmetric  Taylor  vortex  mode  of  instability. 
At  infinite  values  of  Re  solid  body  rotation  is,  according  to  Rayleigh’s  criterion,  neutrally 
stable  and  this  is  the  origin  of  the  real  values  of  K  which  that  eigenvalue  problem  has  in 
the  limit  Re  — »  oo.  We  restrict  our  attention  to  the  case  when  the  pivot  point  for  the 
cylinder  is  midway  along  its  axis.  From  Gerber  et  al.^  we  know  that  the  forced  problem 
at  small  amplitudes  then  leads  to  the  following  conditions  for  u,  v  and  w  at  the  walls: 

u  =  -t(l  -  f)x, 

v  =  -(l-f)x,  (1.9a,  b,  c) 

ta  =  t(l  —  f)r. 


Following  HSG  we  transfer  the  inhomogeneous  boundary  conditions  to  the  endwall  by 
writing 


u  =  — *[1  -  /] i  +  u(r,  x), 
v  =  -[1  -  f]x  +  u(r,x), 
w  =  i[l  -  f]r  +  4>(r)  +  w(r,x), 
P=  ~[1  -  f2}™ +  p(r,x), 


(1.10a,  6,  c,  d,  e) 


4>(r)  =  2  if 


r 


h (MY 

MX) . ' 
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Here  J\  denotes  the  order  1  Bessel  function  of  the  first  type  while  A  is  defined  by 


X  =  f)Re/2. 


(1.11) 


Since  we  will  be  interested  in  the  limit  of  high  Reynolds  numbers  we  can  drop  the  Bessel 
function  terms  in  (1.10)  away  from  r  =  1.  We  note  in  passing  here  that  A,  and  hence  u 
etc.  have  a  branch  cut  at  /  =  1;  we  will  not  investigate  such  high  frequencies  here  since 
the  ballistic  applications  of  our  work  typically  correspond  to  values  of  /  less  than  .5  .  In  a 
separate  study  eigenvalues  and  the  pressure  coefficient  cp  were  calculated  for  1  <  /  <  1.5. 
The  results  were  to  be  compared  with  experimental  data  for  /  >  1.  In  terms  of  u,  v,  w  and 
p  it  follows  from  above  that  the  appropriate  boundary  conditions  are 

u  —  iv  =  w  =  p  =  0,  r  =  0 

u  =  v  =  w  =  0,  r  =  1  (1.12a,  b,  c) 

u  =  v  =  0,w  =  —<f>(r)  x  —  ±A. 


Furthermore,  we  note  that  u,v,w,p  satisfy  (1.5)  and  (1.6)  and  following  the  spatial  eigen¬ 
value  approach  of  HSG  we  therefore  write 


n=l 


q„  sin  knx 
sin  knA 


«n(r). 


Qn  cos  knx 
sin  knA 


Wn(r), 


•-£ 

n=  1 


an  sin  knx 
sin  knA 


”n(r), 


n=  1 


a„  sin  k„x 
sin  knA 


Pn(r). 


(1.13a,  b,  c,  d) 


Here  {&„}  is  the  infinite  sequence  of  eigenvalues  of  (1.6),  (1.7)  and  the  coefficients  {an} 
are  to  be  determined.  In  HSG  these  constants  were  found  by  collocation  at  the  endwalls; 
here  we  shall  make  the  further  limit  Re  — ♦  oc  and  show  how  asymptotic  expressions  for 
these  coefficients  can  be  found. 


2.  The  spatial  eigenvalue  approach  in  the  limit  Re  — ►  oo. 

On  the  basis  of  our  experience  with  the  eigenvalue  problem  for  {fcn}  at  finite  values 
of  Re  we  know  that  the  eigenvalues  split  into  groups  of  three  which  can  be  ordered  by, 
for  example,  connecting  them  to  the  eigenvalues  of  a  related  Taylor  vortex  problem.  We 
anticipate  this  split  and  therefore  consider  the  three  sets  of  eigenvalues  {kjn  },  for  j  —  1, 2, 3 
and  n  =  1,2,3 .. ..  The  distribution  of  eigenvalues  is  illustrated  in  Figure  1  for  the  case 
Re  =  1000  and  /  =  .1;  this  is  taken  from  HSG  wherein  the  labeling  of  the  eigenvalues  is 
discussed.  (Note  that  the  eigenvalue  £3,1  was  labeled  on  the  basis  of  an  analysis  pertinent 
to  the  limit  Re  — ►  0.  However,  in  the  limit  Re  — >  00,  it  lies  in  the  set  of  values  determined 
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by  Stewartson’s  analysis^.)  In  Figure  2  we  show  the  development  of  the  spectrum  as  Re 
increases,  for  fixed  /  and  j  =  1.  We  note  that  for  sufficiently  small  n  the  eigenvalues 
approach  the  real  axis  while  for  large  n  the  eigenvalues  approach  the  imaginary  axis. 

On  the  other  hand  numerical  calculations  in  HSG  showed  that  the  j  =  2,3  branches 
had  |  kjn  j~  Re  *  for  Re  »  1.  We  shall  now  describe  the  precise  details  of  the  asymptotic 
structure  of  {£,„}  for  j  =  1, 2, 3  in  the  limit  Re  -*  oc  and  hence  put  ourselves  into  a  position 
to  find  asymptotic  forms  for  {an}.  Splitting  of  the  eigenvalues  into  triplets  leads  naturally 
to  the  definitions  of  branches  1  ,  2  and  3  in  the  complex  k  plane  for  j  =  1, 2, 3  respectively. 
The  corresponding  amplitude  functions  are  denoted  by  en,7„  and  f3„  ,  replacing  the  an. 

First  we  derive  the  asymptotic  structure  of  the  j  =  1  modes;  the  analysis  is,  of  couise, 
essentially  inviscid  and  based  on  that  given  by  Stewartson^.  We  seek  an  asymptotic 
solution  for  the  nth  eigenvalue  kin  by  writing 

kin  =  k°ln  +  k\nRe~*  +  klnRe~l  +  . . . ,  (2.1) 

while 

«m  =  tt?„  +  u\nRe~  *  +  u\nRe~ 1  +  . . . ,  (2.2) 

together  with  similar  expansions  for  uin,^in  and  pi„.  If  we  define  a  by 


a  =  i{f  -  1} 


then  substitution  of  (2.1),  (2.2)  into  (1.6)  with  I\  replaced  by  kin 
equations  for  u\n,v\n,w\n,p\n  which  has  the  solution 

p2 + 4  k„ = .{(l  -  nA.  +2Jl(f"r)}. 

p2  +  4K„  =  {(1  -  f)M&n  +  2MUM„r)J, 

-k° 

w°ln  =  -f-MAn  r). 

Pin  =  MP°inr). 


yields  a  system  of 


(2.3a,  6,  c,  d) 


/  0  \2  _  r  i.O  \2  (3  “/)(!+  /) 

\Pln)  ~  K^ln)  ^  _  jy 

The  radial  velocity  component  above  satisfies  the  sidewall  condition  if 


(2.4) 


(l-f)PlnJ'Mn)  +  2JMn)  =  0 


(2.5) 
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which  has  an  infinite  number  of  eigenvalues  on  the  positive  real  axis.  [Information  on  this 
spectrum  can  be  found  in  ‘Conduction  of  Heat  in  Solids’  by  H.S.Carslaw  and  J.C.  Jaeger, 
Oxford  University  Press,  1947.]  The  eigenvalues  with  j  =  1  become  progressively  closer 
to  the  real  axis  in  the  limit  Re  — ►  oo.  The  Stewaxtson  eigenvalues  correspond  to  a  choice 
of  frequency  such  that  k°nA  =  y{l  4-  Ar}  for  N  an  even  non-negative  integer.  Thus,  at 
the  resonant  frequencies  of  a  cylinder  of  given  length,  the  j  =  1  eigenvalues  approach 
Stewartson’s  eigenvalues.  Since  we  have  not  retained  viscous  effects  in  the  perturbation 
equations,  we  can,  of  course,  only  make  the  radial  velocity  component  vanish  at  r  =  1. 
The  other  velocity  components  must  therefore  be  reduced  to  zero  within  a  boundary  layer 
of  thickness  Re~t  near  r  =  1.  The  determination  of  the  velocity  field  within  this  layer 
will  enable  us  to  determine  the  order  Re~3  term  in  (2.1),  we  shall  give  some  of  the  detail 
necessary  to  determine  it. 

First  we  note  that  the  perturbation  equations  for  r  =  0(1)  proceed  in  powers  of  Re  so 
that  if  we  define  /iln  by 

Pin  =  Pin  +  Pin  Re  2H - 5  (2-6) 

then  correct  to  order  Re~  2,  (2.3)  are  valid  if  u°n1v^n,w°n,p°n  and  p°n  are  replaced  by 
uin  +  uin  etc.  and  then 


„0  1  _kQ  (3  -  /X*  +  /) 

P lnPln  ~  ^  f-y 


It  follows  that  when  r 


(2.7) 


(3  -  /)(1  +  f)u\n 


~Pln  “ /)  1  + 


(3  -  /)(!  +  /) 

K„)2(W)2 


(2.8) 


where  (  =  Re  3  (1  —  r),  and  C  is  a  constant  to  be  determined  later.  In  the  boundary  layer 
at  r  =  1  we  write 


Re3Uln  =  u°ln  -f  u\nRe  2  H - 

together  with  similar  expansions  for  ui„,u;in,  and  pJn.  It  can  be  shown  that  the  zeroth 
order  system  obtained  by  substituting  these  expansions  into  (1.6)  written  in  terms  of  £ 
yields 

=  constant,  ^  }  {l  +  (*!„)2}  {  ^  ^  ^  +  C> 
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when  solved  subject  to  =  v°n  =  w°n  =  0,  £  =  0.  Thus  matching  the  pressure  and 
radial  velocity  component  when  the  inviscid  and  viscous  regions  meet  give 

-o  _  r  ,  o  ,  iJ.(A)|l +  (*?„)*]  _  -  /W,  r,  ,  (3 -/)(!+/)] 

Pl"  1.).  (1_/Vi  (3-/)(l+/)  t  +«)2(l-/p/’ 

and  an  expression  for  C  the  constant  in  (2.8).  Using  (2.7)  and  the  eigenrelation  (2.5)  the 
second  of  the  above  equations  can  be  manipulated  to  give 

^  =  +  +  (2'9) 


where  k°n  is  determined  by  (2.4)  and  (2.5). 

The  viscous  modes  have  an  x-dependence  on  a  Re~  2  lengthscale.  For  convenience  we 
drop  the  ‘2n’  and  ‘3n’  notation  for  the  moment  and  seek  a  viscous  eigenvalue  of  (1.6)  and 
(1.7)  by  writing 

k  =  y/Re  ko  H — 7=^1  H -  (2.10) 

y/Re 

while  u,  v,  w,  and  p  expand  as 


u  =  Uq  + 


\/Re 


u  1  + 


v  =  Vo  H - rvi  H - 


w  = 


Wq 


Re i 

W\ 

+  ^-  + 


y/Re  i?e 


„  _  Po  Pi 

The  zeroth  order  approximation  to  (1.6)  then  yields 

{— ‘kg  —  if  +  i}u0  +  2i>o  =  0, 

{-&o  -  if  +  i}v 0  -  2u0  =  0, 


du0 

W°  =  -d//k°- 


The  equations  (2.12a,b)  are  consistent  if 


kl  =  (3  -  /).' 
or  k20  =  -(1  +  j)i 


(2.11a,  b,  c,  d ) 


(2.12a,  6,  c) 


(2.13a,  6) 
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which  correspond  to  the  j  =  3  and  j  =  2  modes  respectively.  At  higher  order  uo>  and  hence 
t>o,u>o  and  Poi  are  determined  as  solvability  conditions.  We  find  that  (2.13a,b)  respectively 
lead  to 

u0  =  J2(snr),  u0  =  Jo(tnr)  (2.14a,  6) 


where 


Msn)  =  J0(tn )  =  0. 


A  more  detailed  derivation  of  (2.14)  is  found  in  HSG.  Thus  (2.13a,b)  correspond  to  two 
countably  infinite  sets  of  eigenvalues  and  {<sn},{tn}  determine  the  order  Re~?  correction 
terms  in  the  expansion  of  k.  In  fact  the  appropriate  equations  axe  found  to  be 


and 


h  = 


*i  = 


-(1-0(4-/)^ 
2^2(3  -  7f 

~(2 +  /)<„(!  + 1) 


(2.15) 


(2.16) 


2^2(1  + /)3/2 

Henceforth  we  refer  to  the  eigenvalues  corresponding  to  (2.13a,b)  using  the  indices  j  =  3 
and  j  —  2  respectively.  The  structure  outlined  above  breaks  down  when  n  is  formally 
0 (Re)  in  which  case  the  eigenfunction  has  both  axial  and  radial  variations  on  the  Re~* 
lengthscale.  We  can  now  modify  the  spatial  eigenvalue  expansions  (1.13a,b,c,d)  using  the 
asymptotic  structure  developed  above;  we  therefore  write: 


eni 


(3  —  /}{!+/} 


{ 


(/  +  l)*(#*?„r) 


+  (1  -  /)m 


o  t  /  o  a!  sinfc?n 

-■7o('‘-r)/^ 


x 

A 


+  7  nJo(tnr)~—j  +  /3nJ2(snr)~Y~2 

smAiA 


v=E 

l  L 


sin  A  2  A 
Cn 


{3 -/}{!+/} 


a  ,  r)  ,  „  0  r  / .  o  sin k°nx 

"(1  +  f)^—  +  2^J^r)  | 


r  sin  A2x  ,  T  ,  ,  sinAjx 

-  i-ynJo{tnr)-. — -  7  4-  ipnJ2{snr)  — 


w 


sin  A2A 

E°°  .  .  o  .COsfcfng 

;  (/ _ 


sin  Ai  A 


+ 


p  =  Y^enJMnr)-: 

*  Si 


sin  fcJnA 

sinfcjnx 
sin  fc°„A 


(2.17a,  6,  c,d) 
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where 


Aj  =(1  +  0 


3-/ 


y/Re,  A2  =  (1  —  0 


1  +  / 


\/Re. 


We  now  apply  the  endwall  condition  (1.12c)  on  the  axial  velocity  component  w  to  give 


—2  ifr  =  £ 


kl»A 


so  that  after  multiplying  both  sides  of  the  equation  by  rj\(p({nr),  integrating  from  r  =  0 
to  r  =  1  and  using  the  well-known  properties  of  Bessel  functions  we  obtain 


— 


4/(3  -  /)  tan  k'fnA 

+<±?)Mrfn) 


(2.18) 


while  the  endwall  conditions  on  the  radial  and  circumferential  velocity  components  jield 
equations  which  can  be  solved  for  {/?„}  and  {7„}  in  terms  of  {en}.  It  is  at  this  stage 
that  the  differences  between  the  finite  and  infinite  Reynolds  number  have  emerged.  Thus 
the  amplitude  functions  {e„},  {/3n}  and  {7,,}  in  the  limit  Re  approaches  00  can  be  found 
explicitly  using  orthogonality  properties  of  Bessel  functions.  At  finite  values  of  Re  HSG 
found  it  necessary  to  use  either  collocation  or  a  least  squares  approach  in  order  to  find  the 
corresponding  amplitudes.  At  higher  order  in  the  procedure  that  we  have  begun  above, 
the  sidewall  boundary  layer  comes  into  play  at  the  endwalls;  at  that  stage  the  higher  order 
amplitude  coefficients  must  be  found  by  collocation  or  the  least  squares  method. 

It  follows  from  (2.17d)  and  (2.18)  that  p  on  the  endwall,  z  =  A,  can  be  expressed  as 
the  infinite  series 


P  =  4/(3  —  /)£ 


1 


tan  k°nA  Ji(p°lnr) 


(2.19) 


which  is  correct  up  to  order  Re0. 

The  pressure  coefficient 

We  shall  now  use  the  high  Reynolds  number  form  of  the  spatial  eigenvalue  approach 
to  find  the  pressure  exerted  by  the  fluid  on  the  container.  Since  almost  all  of  the  available 
experimental  results  for  pressure  correspond  to  measurements  made  at  the  endwall  we  shall 
concentrate  on  that  situation.  More  precisely  we  use  our  theory  to  find  the  amplitude  of 
the  pressure  measured  by  a  pressure  transducer  located  at  the  endwall.  More  detail  of  the 
derivation  of  this  pressure  coefficient  can  be  found  in  Gerber  et  al.W 

First  we  note  that  p+  ,  the  pressure  calculated  from  the  Navier-Stokes  equations  is 


P+  =  \r2  -  KoP* , 
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in  the  case  Kq  — »  0.  We  must  now  relate  our  calculations  in  an  inertial  frame  to  the  frame 
r,9,x  used  in  HSG  where  the  x  axis  coincides  with  the  cylinder  axis.  A  routine  calculation 
based  on  a  Taylor  series  expansion  shows  that,  for  example 


\r2  ~  “  K°fi  cos{-f*  -&}  +  0(A'o ). 

We  now  let  A p  be  the  disturbance  pressure  not  including  the  contribution  from  solid  body 
rotation.  At  points  fixed  on  the  surface  of  the  cylinder 

A p=^r2  -  Kop*(r,0,x,t)  -  ^r2 

and  since  p*(r,6,x,t )  —  p*(r,0,x,t)  =  O(ATo)  it  follows  that 

A  p  =  -Kol-jy  sin  (ft  -  0)  +  ( pr  +  rx)  cos  (ft  -  9)}  +  0(A'q  ). 


Here  subscripts  r  and  i  denote  the  real  and  imaginary  parts  of  a  complex  quantity.  Now 
since 

pr  +  rx  =  pr  +  frx,  p.  =  pi 

the  pressure  coefficient  cp  =  may  be  written  as 

K0 

cp  =  ^  (pr  +  rx)2  +p].  (2.20) 

It  then  follows  from  (2.19)  and  (2.20)  that  cv  calculated  using  the  spatial  eigenvalue  ap¬ 
proach  in  the  limit  Re  — »  oo  is  given  by 


cp  — 


4/(1  -  / y 


i +/ 


E 


tan  k%„A  JM„r) 

1  + 


+  fr  A 


+  0(Be-i). 


(2.21) 


The  above  equation  is  not  uniformly  valid  since  at  the  Stewaxtson  eigenfrequencies 
the  unperturbed  state  can  support  neutral  eigenfunctions.  However,  it  is  clear  that  cp 
will  become  unbounded  if  kyn  corresponds  to  one  of  the  Stewartson  eigenvalues  in  which 
case  tan  k°nA  =  oo.  In  fact  our  expansion  procedure  is  readily  modified  to  take  account 
of  this  possibility.  Suppose  then  that  /  is  an  eigenfrequency  f  of  the  inviscid  problem. 
Without  any  loss  of  generality  we  can  suppose  that  fcjjA  =  tt/2  so  that  is  formally 
infinite.  In  order  to  remove  this  singularity  we  note  from  (2.18)  that  if  higher  order  terms 
are  retained  then  ci  is  in  fact  O(-Res)  when  /  =  /.  Since  k\n  =  k\n(f )  and  we  wish  to 
obtain  a  correct  zeroth  order  approximation  to  cp  for  all  values  of  /  we  must  therefore 
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modify  (2.17)  for  /  —  /  ~  0(J?e-  2  ).  Thus  (2.21)  is  the  correct  zeroth  order  approximation 
to  cp  for  f  —  f  >>  0(fZe~2);  we  refer  to  this  situation  as  the  nonresonant  case. 

It  follows  from  our  calculations  above  for  the  nonresonant  case  by  talcing  the  limit 
f  —*  f  that  ei,  {/?„},  {7„}  for  all  n  becomes  formally  0(fJe+2)  when  /  — /=  0(iJe-2). 

We  therefore  make  use  of  the  new  scalings  inferred  above  and  therefore  replace  (2.17) 

by 


ieiVRe  f(l  + 
ti“(3-/)(l+/)l  r 

oo 


+  Z!  (3-/)(l  +/)(  r!  Jl  +  (1  /)^"J°(^r)  +  "}sinfc0^1 

+  VRep,MUr)^  +  - 

+  *U±0.M-r)gfa+... 

eiy/Re  f  „  o  T/  o  nn  .  1  sin kux 

v  =  (3-/K1  +  /){-(1  +  })— i—  +  2t,“' 7«(""r))  +  •  ■ 

+ 1  *■*«•'> + ' 

W  =  +  -]^f 

I  V  ikln  _  r  t  /  0  i  ,C06fcJnI 

+  y-'  snPnj'2(snr)(l  -  i)  cosAix  ^ 

“  \/3  -  fy/ 2  sin  Aj  A  + 

+  y  -7ntnJi(tnr)(l  -ft)  cosA2x  ^ 
j  >/l  4-  /\/5  sin  A2A 

P  =  £lvKj1(/i;,r)  |2ii2£+0(i), 

(2.22a,  6, 


•  +  ••• 


+ 


11 


We  note  that  in  the  above  expansion  /  formally  differs  from  /  by  0 (Re~*).  If  we 
now  apply  the  endwall  conditions  on  the  radial  and  circumferential  velocity  components 
we  find,  after  using  the  well-known  orthogonality  properties  of  Bessel  functions,  that 


W?((«)  =  l'  rMAr)Mt.r)dr, 

The  endwall  condition  (1.12c)  on  the  axial  velocity  component  then  gives 

J1(/j,°llr)cot  kuA  +  cot 

Sn0 it^2(^t»T’)(^  0  C°^  ^lA 

+  2^  ./o./TZT 


(2.23a,  b). 


5n/^n‘^2(5r»7’)(^  i)  COt  \\A 
r-N  ~7n^n  «A(fnr )(1  ~f  Q  COt  A2A 

+  ^  V2VTT7 


=  —2  */r. 


We  now  multiply  the  above  equation  by  rJi(/ijjr)  and  integrate  from  r  =  0  to  r  =  1 
and  replace  /?„  and  7„  using  (2.23).  After  some  manipulations  we  can  show  that 


VnTe 

2 


(1  — l)cotA2A  „  2  j2f  0  ] 

■  V5(1  +  /)’/»  ("n)  •7”('‘“)52J 

-2i/  .  o  X 

~  ,,0  ^2(^11) 

Mil 


(2.24) 


where 


~  _  Vp  _ ■Sn^n _ 

oo  ^2 

Si  = 

Here  zn  =  rj'2(snr)Ji(fJ,iir)dr.  We  note  that  in  (2.24)  y/RecotkuA  ~  0(i2e°)  because 
/  differs  from  /  by  0(Re~  * ).  Finally  we  note  that  cot  Ai  A  and  cot  A2A  to  this  order  can  be 
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replaced  by  —i  and  +i  respectively.  Thus  with  ei  given  by  (2.24)  the  pressure  perturbation 
is 

p=e1y/teJ1(»°nr)  (2.25) 

which  is  of  course  0(ile2)  larger  than  that  in  the  non-resonant  case.  It  can  be  shown  from 
(2.24)  that  when  1  >>/  — /  >>  Re~i  the  equation  (2.25)  is  consistent  with  that  obtained 
from  (2.19)  by  talcing  the  limit  /  — ►  /.  Thus  the  formula 


.  v/n 7  7  /..o  ,  ^4/(l-/)2  tan  k\nA  J\(p\nr) 


(2.26) 


gives  a  correct  zeroth  order  approximation  to  p  for  all  /.  The  pressure  coefficient  cp  is 
thus  given  by 


cp  ~ 


4/(1  -  /)2 


E 


tan  k°lnA  Ji(p°lnr) 


1  +  /  V 

+f2rA  +  ei\/^eJi(/x?1r)| 


(2.27) 


with  €i  given  by  (2.24). 

We  close  this  section  with  a  discussion  of  the  relationship  between  our  resonant  expan¬ 
sion  procedure  discussed  above  with  the  viscous  correction  to  Stewartson’s  inviscid  theory 
given  by  Wedemeyer^ .  In  the  latter  paper  Wedemeyer  showed  how  a  viscous  correction 
to  the  inviscid  theory  could  be  obtained  by  considering  viscous  boundary  layers  at  the 
ends  and  sidewall  of  the  cylinder.  (We  observe  that  a  and  (3  of  Eqs.(33a)  and  (33b)  of 
Reference  6  are  related  to  Aj  and  Aj  on  p.9  of  this  paper  as  follows:  Ai  =  iaa,  A2  =  —ia)3.) 
We  note  that  away  from  the  endwalls  the  terms  that  appear  inside  the  summation  signs  in 
(2.22)  are  exponentially  small.  The  term  that  remains  in  the  expressions  for  the  velocity 
and  pressure  corresponds  to  Stewartson’s  eigenvalue.  However  the  axial  wavenumber  used 
there  is  &u  correct  to  order  Re  1/^2.  Thus  there  is  a  viscous  correction  in  our  expansion 
that  arises  from  the  fact  that  we  choose  to  continue  the  inviscid  expansion  in  order  to 
take  care  of  the  viscous  sidewall  layer.  Thus  the  modified  eigenfunction  we  have  used  is 
precisely  that  implied  by  Wedemeyer’s  equation  (16).  At  the  ends  of  the  cylinder  all  the 
ratios  of  the  trigonometric  functions  are  0(1)  so  that  all  the  terms  are  to  be  retained.  Now 
the  terms  proportional  to  sin(Xjx),  cos(Xjx)  are  needed  in  order  to  satisfy  the  endwall 
conditions.  If  we  choose  to  look  at  the  unforced  problem,  then  it  can  be  shown  that  the 
effect  of  these  terms  is  exactly  equivalent,  in  a  formal  asymptotic  sense,  to  Wedemeyer’s 
conditions  (18,19).  Thus,  for  the  unforced  problem  at  the  resonant  frequency,  our  expan¬ 
sion  is  equivalent  to  that  of  Wedemeyer.  Notice,  however,  that  the  form  of  our  expansion 


13 


enables  us  to  determine  in  a  precise  asymptotic  form  the  forced  velocity  field  both  at  the 
resonant  frequency  _nd  elsewhere.  That  is  not  the  case  for  Wedemeyer’s  theory  so  we  can 
interpret  (2.22)  as  being  a  generalization  of  Wedemeyer’s  work  to  the  forced  problem  at 
all  frequencies. 

3.  Results  and  Discussion 

We  shall  give  results  only  for  the  composite  expansion  (2.27)  rather  than  plotting 
separately  the  resonant  and  non-resonant  expansions.  In  particular  we  will  compare  our 
results  with  those  of  Whiting^  who  investigated  the  problem  experimentally.  In  Figures 
3a,b,c  we  compare  with  experimental  data  for  the  pressure  coefficient  measured  at  r  =  .668 
with  A  =  3.1481  at  three  different  values  of  the  Reynolds  number.  Also  shown  in  these 
figures  are  curves  representing  theoretical  results  from  Gerber  et  al.M  and  Murphy^.  Both 
of  the  latter  calculations  used  approximate  endwall  conditions  obtained  using  Wedemeyer’s 
approach.  At  Re  =  5  x  103  the  results  of  the  present  approach  and  Reference  7  agree  more 
closely  with  experiment  than  the  results  of  Reference  4  theory.  At  Re  =  105  the  present 
theory  gives  as  good  agreement  with  the  theory  as  the  other  approaches  but  at  the  highest 
value  of  the  Reynolds  number  our  method  overestimates  the  pressure  coefficient  at  the 
resonant  frequency.  Apart  from  the  immediate  neighborhood  of  the  resonant  frequency 
the  different  theoretical  approaches  at  the  two  highest  values  of  Re  are  in  agreement. 
The  difference  between  the  predicted  resonant  response  and  the  experimental  data  could 
presumably  be  eliminated  by  proceeding  with  the  present  calculation  to  the  next  order. 

In  Figure  4  we  show  a  similar  comparison  at  A  =  1.0509;  again  we  see  that  the 
differences  in  the  two  approaches  are  small  except  for  a  small  interval  near  the  resonant 
frequency.  In  Figure  5  we  show  the  variation  of  the  resonant  frequency  with  Reynolds 
number  for  A  =  3.1481. 

In  Figure  6  we  show  3D  plots  of  the  pressure  coefficient  at  two  different  Reynolds 
numbers  over  a  portion  of  the  aspect  ratio-frequency  plane.  Equal  intervals  in  /  and  A  are 
used  in  these  plots.  The  localized  peaks  correspond  to  resonant  responses  caused  by  the 
excitation  of  one  of  the  Stewartson  modes.  Such  plots  are  useful  in  picking  out  portions  of 
parameter  space  with  certain  properties;  more  detailed  plots  or  variations  can  be  obtained 
from  the  program  written  to  calculate  cp.  A  few  points  about  the  approach  that  we  have 
developed  here  are  in  order: 

1.  The  composite  expansion  is  relevant  whenever  any  one  of  the  Stewartson  eigenvalues 
occurs;  the  eigenvalue  whose  eigenfunction  is  rescaled  by  a  factor  \/Re  is  taken  to  be  that 
eigenvalue. 
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2.  The  computational  time  for  the  present  approach  is  comparable  to  that  required 
for  the  different  approaches  mentioned  above.  For  example  all  the  calculations  used  to 
supply  the  data  for  the  3D  plots  (with  2500  nodes)  were  carried  out  in  about  one  hour  on 
an  Olivetti  M24  PC. 
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1.  The  first  72  eigenvalues  for  Re  ■  1000,f  «=  0.1 
(taken  from  Reference  2). _ 
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Figure  2.  The  eigenvalue  spectrum  for  some  large  values  of 
and  f  *  0.08 


Figure  3a.  The  behavior  of  Cp  calculated  using  the  composite 

expansion  with  A  ■  3.1481,  Re  =  5000.  The  experimental 
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Figure  4a.  The  behavior  of  cp  calculated  using  the  composite 
expansion  with  A  *  1.0509,  Re  =  100,000. _ 
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Figure  4b.  The  behavior  of  cp  calculated  using  the  composite 
expansion  with  A  =  1.0509,  Re  -  500>000. _ 
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Figure  5.  Variation  of  the  frequency  of  peak  response  with  Re 
for  A  =  3.148.  Solid  line  is  from  the  approach  of 
Reference  4.  The  crosses  are  from  the  present  cal¬ 
culations  . 
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Figure  6b.  A  3D  plot  of  Cp  against  A  and  f  for  the  case  Re  =  500,000. 


Intentionally  left  biank. 
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